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i , Abstract 

We discuss possible phase factors for the S-matrix of planar Af = 4 gauge 
theory, leading to modifications at four-loop order as compared to an earlier 
proposal. While these result in a four-loop breakdown of perturbative BMN- 
scaling, Kotikov-Lipatov transcendentality in the universal scaling function 
for large-spin twist operators may be preserved. One particularly natural 
choice, unique up to one constant, modifies the overall contribution of all 
terms containing odd zeta functions in the earlier proposed scaling function 
based on a trivial phase. Excitingly, we present evidence that this choice 
is non-perturbatively related to a recently conjectured crossing-symmetric 
phase factor for perturbative string theory on AdS^ x S 5 once the constant 
is fixed to a particular value. Our proposal, if true, might therefore resolve 
the long-standing AdS/CFT discrepancies between gauge and string theory. 



1 Introduction 



Supersymmetric M = 4 gauge theory is a conformal quantum field theory in four dimen- 
sions. As such, standard lore says that it does not possess an S-matrix since there are 
no asymptotic particles. Luckily, this naive "no-go" theorem is wrong in two interesting 
ways. 

Firstly, it turns out that an asymptotic "world-sheet" S-matrix of this theory may 
be defined, in the planar limit, in an internal space [T]. This is possible since planar 
local composite operators may be interpreted as a one-dimensional ring on which 8 
elementary bosons and 8 elementary fermions propagate [2]. Furthermore, this S-matrix 
appears to be two-particle factorized, since the particle model is integrable, at one-loop 
and beyond [3]. This allows, under some assumptions, to largely construct the full 
S-matrix |4J, up to an a priori unknown phase factor, and find the spectrum of the 
theory [5]. This phase factor is known to equal one up to three-loop order. 

Secondly, it is very interesting to consider the set of all space-time on-shell n-gluon 
amplitudes of the Af = 4 gauge theory. While these are clearly not scattering amplitudes 
of asymptotic particles, they contain much physical information, and we can simply 
declare them to serve as a surrogate "space-time" S-matrix. Fascinatingly, much evidence 
for an iterative "solvable" structure in these amplitudes was discovered recently [HIE]- 

The factorized world-sheet S-matrix and the iterative space-time S-matrix approaches 
are not uncorrelated. The conjectured ansatz for the all-loop gluon amplitudes contains 
a to-be-determined function f(g), where 



is the gauge coupling constant. However, this function f(g) may be extracted from the 
logarithmically divergent anomalous dimension of leading-twist operators in the gauge 
theory [8]. As such, it may in principle be derived from the diagonalization of the world- 
sheet S-matrix. In this work we will call f(g) the universal scaling function of M = 4 
gauge theory; it is also known as the cusp or soft anomalous dimension. 

In practice, the comparison has so far been done reliably to three-loop order only. 
The space-time based approach was completed in [7j, and the world-sheet method was 
exploited in pp. The results agreed. The latter method was also successfully checked 
against rigorously known anomalous dimensions of twist operators at two loops, and 
brilliantly conjectured ones at three loops [9J. The conjecture was based on an intrigu- 
ing transcendentality principle, which allowed to extract the answer from a hard QCD 
computation [TO] . 

Clearly, one would like to go to higher loops, and ideally compare the full univer- 
sal scaling function. A four-loop computation by Bern, Czakon, Dixon, Kosower and 
Smirnov using the iterative gluon amplitude approach is under way [Hjo Furthermore, 
in [13] it was shown how to extract the universal scaling function from the Bethe equa- 
tions which result from the diagonalized S-matrix. Assuming a trivial (i.e. equal to unity) 
phase factor, an all-loop equation determining the function was proposed. While the ex- 
act solution of the equation remains unknown, it allows to generate perturbative results 

1 While finalizing our manuscript we were informed that this computation [12] has been completed. 
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to any desired order with ease. In addition, the solution preserves the Kotikov-Lipatov 
transcendentality principle [14J to arbitrary loop order. 

However, one troubling issue is whether the phase factor is really unity to all orders in 
perturbation theory. This leads to various vexing discrepancies [ToUTO] with string theory 
results in the context of the AdS/CFT correspondence [17]. In fact, the S-matrix [4"l[18] 
of string theory on AdS^ x S 5 definitely requires a phase factor pj5] as can be seen 
by comparison to an integral equation describing classical solutions [20]. Furthermore, 
it was argued that the S-matrix should possess a crossing symmetry [21], which then 
necessarily calls for such a phase. In [22] it was shown that the known classical and one- 
loop part of the phase factor on the string side satisfies crossing symmetry, and in [23] a 
crossing-invariant all-order string phase factor was proposed. 

In the weakly coupled gauge theory, the dressing phase could appear, at the earliest, 
at four- loop order [21] . The leading (four- loop) effect on the scaling function was already 
investigated in [T3]. Here we will extend this analysis to all- loop order, and investigate 
the effects of the dressing phase on the universal scaling function. 

2 Integral Equations for the Scaling Function 

Let us start in this section by considering the effect of the dressing phase on the universal 
scaling function. We first briefly recall the assumptions and main steps leading to an 
integral equation for the scaling function of leading twist operators of M = 4 gauge 
theory in the large spin limit S — > oo. Up to three loops, it is based on the asymptotic 
Bethe ansatz 



This ansatz was shown to properly work in the sl(2) sector relevant to twist operators 
up to three- loop order in [I], and further tests were performed in [25J. It was constructed 
to possibly describe gauge theory at four loops and beyond in |5], in analogy with an 
earlier proposal for the su(2) sector [26] . 

The large-spin computation is similar to a thermodynamic limit for a spin chain 
Bethe ansatz, as the number S of Bethe roots gets very large, and their distribution may 
be described by a smooth density. Note, however, that the actual length L of this spin 
chain should remain quite short. Ideally it should be L = 2, but this is dangerous since 
the Bethe ansatz is a priori only asymptotic [HUES, 1, 5J. A work-around was devised 
in [13], where it was argued that the scaling of the anomalous dimension of the lowest 
state of operators of finite twist-L is a universal function of the coupling constant q. This 
suggests that it should coincide with the scaling function of twist-two operators^ For a 
similar discussion, albeit at one loop, see [27J. 

2 At four loops the asymptotic Bethe ansatze are only fully reliable at twist three, or higher, and we 
have to assume the four-loop universality of the scaling function for the lowest state of twist operators. 
This universality may only be rigorously proved for twist > 2 at one loop [27] , for twist > 2 at two and 
three loops, and for twist > 3 at four loops [T3] . It would clearly be most welcome to more directly 
confirm universality from the field theory. 
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Despite the similarity to a thermodynamic limit, the computation is quite subtle due 
to the singular distribution of the one-loop roots [28] . The problem was solved in [13] by 
splitting off the one-loop piece, and subsequently deriving an integral equations for the 
non- singular fluctuations around it. It was found that the density of fluctuations a(t) is 
determined by the solution of a non-singular integral equation of the general form 

a(t) = -— ■ K{2 gt,0)-4g 2 dt' K(2 gt,2g f) a{t') . (2.2) 
e — i. i J J 

The universal scaling function f(g) is then given by 

f{9) = 8g 2 -64g* ^ dta{t) J 4r^- (2-3) 
Jo 2 9t 

Here and in the following J r (t) denotes a standard Bessel function. Note that this 
expression can be reduced provided that the kernel satisfies the property (it will always 
be true in this article) 

K(0,t) = ^. (2.4) 

The scaling function then simply takes the value of the fluctuation density at t — [29] 

f(g) = 16 g 2 &(0). (2.5) 

For the above, particular Bethe ansatz ( 12 .ip the kernel K(t,t') = K m (t,t') is given 
by (m stands for main scattering) 

(. ,A -hit) ut') - ut) .hit') 

The inhomogeneous piece of the integral equation ( 12.21) with this kernel thus reads 

KUt,0) = ^. (2.7) 

For this specific kernel (12. 6p derived from (12. ip the weak-coupling expansion of the scaling 
function reads 

U9) = *9 2 -l *V + § vr V " 16 vr 6 - 4 C(3) 2 ) / 

+ 32 (J?L n s _ ^ ^(3)2 _ 4Q c(3) c(5) ^j 5 io 

V3742200 15 SV ; 3 SV ;sv ; 

-210C(3)C(7)-102C(5) 2 ^) g 12 

G 



/ 7680089 12 _ _4J_ 6 2 _ 82 4 _ ^2^3) c(7) 

1 340540200 189 sv ; 15 sv ;sv ; sv ;sv ; 



-34tt 2 C(5) 2 - 1176C(3)C(9)- 1092C(5)C(7) + 4C(3) 4 J g 14 
T ■ ■ ■ ■ (2.8) 
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Here we can observe the Kotikov-Lipatov transcendentality principle [H]: We attribute 
a degree of transcendentality k to the constants n k as well as to C(&). 

The i-loop contribution to the universal scaling function f(g) forM = 4 gauge theory 
has a uniform degree of transcendentality 2£ — 2. 

The function f m (g) derived from the integral equation (12.21) obeys this rule [13]. A 
rigorous proof to all orders in perturbation theory was given in [29|. 

Let us next generalize the universal scaling function to the case of an arbitrary weak- 
coupling dressing phase. The only restriction is that the phase should merely modify the 
gauge theory Bethe ansatz at four loops or beyond, as the Bethe ansatz (12. ip is firmly 
established up to three-loop order [301124]. The corrections of the higher-loop Bethe 
equations which are consistent with current knowledge are 

\ x k / j = i x k x j 1 y I x k x j 

where the dressing phase 9 is conjectured to be of the general form [HHl3"Tll2"4"] 

oo oo 

6(u k , Uj) = ^2^2 Pr,r+i+2u{y) (q r (uk) q r +l+2u{Uj) - q r {uj) q r +i+2u{u k )) , (2.10) 

r=2 i/=0 

The q r (u) are the eigenvalues of the conserved magnon charges, see [26J. The coefficient 
functions (3 r ,s{y) expand in powers of the coupling constant 

oo 

/WW<7) = Y,9 2r+2 " +2 ^t + + Zi (2-11) 

fX=U 

with /?rfi some numerical constants in the same notation as in [21]. We should eliminate 
from the start the only coefficient f3\ 3 = contributing at three loops. Note that in [13] 
only the leading four-loop correction 

exp(2i6(u k ,Uj)) = exp(i8pg 6 {q 2 {u k ) qz{uj) - q 3 {u k ) q 2 (uj)) + . . .) (2.12) 

was treated. Thus the first constant in fl 2 . X 1) is j3\ 3 = 4/3. 

The modified Bethe ansatz (12.91) with (12.101) may be treated in much the same fashion 
as the simpler ansatz (12.11) . see [13] for the details. In particular, the singular one-loop 
distribution may be split off in the same way. In fact, for an arbitrary dressing phase 
the general form ( 12.21) of the integral equation as well as the properties ( I2.5f2.4l) remain 
valid. What changes is the kernel ( 12.61) . which generalizes to 

K(t,t') = k m (t,t')+k d (t,t') (2.13) 

with the dressing kernel 

. 00 00 00 

p=l u=0 /i=i/ 



-if - WW - Vi(0 (2.i4) 

+ /^2p+1^2p+2^+2 J2p{t) J2p+l+2v(t')) , 
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and the inhomogeneous term of the integral equation ( 12.2ft receives additional contribu- 
tions 

*d(f,o) = ^ i^ErHrC:' -w*) • (2.15) 

Notice that as promised ( 12.41) still holds for ( 12.13112.141) and arbitrary constants (3^1, and 
that therefore also ( 12.51) remains valid. 

Equation (12.21) is still just as suitable for a perturbative small-g expansion if we use 
the modified kernel (I2.13f2.14l) instead of (12. 6p . It was already stated in [T3] that a 
four-loop dressing phase (I2.12p leads to the following 0(g 8 ) term in the scaling function 

/(//): 

M = ■ ■ ■ - 16 Q| vr 6 - 4C(3) 2 + 8/3 C(3)) / + •••■ (2.16) 

This violates transcendentality for generic (3. However, if (3 is a rational number times 
£(3) (or 7r 3 ) transcendentality is preserved. It is straightforward to extend this analysis 
to even higher loop order: 

Interestingly, a study of the effect of the constants (3$ in ( Iff. 1 0\) reveals that for arbi- 
trary loop order transcendentality is preserved if and only if the degree of transcendentality 
of $r+i+2v i s 2/J- + 1, independently of r and v. In other words 

(3^1 should have degree of transcendentality 21 + 2 — r — s. (2-17) 

A particularly curious case is (3 — |C(3), which cancels the term containing C(3) 2 , 
and thus leads to the much simpler four-loop answer — ^ tc 6 16 g 8 . Beyond four loops, 
one finds that one can always choose the 4fs in many different ways such that all terms 
containing zeta functions of odd argument are canceled from the expansion ( 12. 8ft . What 
is truly remarkable, however, is that the constants 4fs are uniquely determined if we 
impose a further restriction on them conjectured to hold for arbitrary long-range spin 
chains compatible with gauge theory^] in [21], cf. eq. (4.2) in that paper, 

= for fJKr + u-l, (2.18) 

or, in a different notation, (3^1 = for £ < r + s — 2. One finds to the first few orders^ 

+2C(3), 
-20C(5), 

+210 C(7), 42 - +12C(5), 45 - -4C(5), 

-2352C(9), 4 6 j^-210C(7), 4$ - +84C(7). (2.19) 

3 This study applied to the su(2) sector of the gauge theory, but the phase is universal for all sectors 
[UG3U and thus includes st(2). 

4 Some of the coefficients appear in the expansion of the scaling function at higher orders than naively 
expected from the expansion of the kernel (|2.13I2. 14|) . They are nevertheless fixed. E.g. f3^\ does not 
appear at six loops as expected, but only at seven loops in f(g). But at this order it is fixed to 12 £(5) 
if we demand all odd-zeta terms to cancel. 
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The expansion of the scaling function significantly simplifies as compared to (12.81) : 

f (g) = 8 g 2 - - vr 2 g* + — vr 4 g 6 - 16 — vr 6 g 8 + 32 n*g 10 
joyyj y 3 y ^ 45 y 630 y 14175 y 

- 64 . V* + 128 7680089 ^ T . . . . (2.20) 

3742200 y 340540200 y v ; 

The even zeta terms, and thus the parts containing only even powers of tt, are unaffected. 

However, (I2.19P is not the only curious choice for the constants /3^} in the dressing 
phase (12.101) . Another striking choice corresponds to doubling the just discussed special 
constants, e.g. to the first few order (12.191) become: 

4 3 3 } = +4C(3), 
/3$ = -40C(5), 

43 = +420C(7), /fj = +24C(5), = -8C(5), 

/?g = -4704C(9), /?g = -420C(7), $j> = +168 C(7). (2.21) 

Now the zeta functions of odd argument no longer cancel out. Instead, one finds to 
e.g. seven-loop order 

M = 8 9 2 ~l vr V + I vr V " 16 ^ + 4 C(3) 2 ) 9* 



+ 32 (J?L ^ + 4 ^2 + 4Q c(3) c(5) ^ gW 

-64fi^^ + AA(3) 2 + ^^ C (3)C(5) 
1 3742200 15 sv ; 3 SV ;sv ; 



+ 210C(3)C(7) + 102C(5) 2 U 



.12 



/ 7680089 12 _47_ 6 2 82 4 70tt 2 C(3) C(7) 

V 340540200 189 sv ; 15 swsw 



+ 34vr 2 C(5) 2 + 1176 C(3) C(9) + 1092 C(5) ((7) + 4 C(3) 4 J ^ 14 

=F • • ■ • (2.22) 

Remarkably, the alternating sum (12.221) is identical to (12. 8ft for the case of a trivial dress- 
ing phase by multiplying all zeta functions with odd arguments by the imaginary unit i, 
i.e. the replacement ((2n + 1) — > i((2n + 1). After this operation, and in contradistinc- 
tion the the earlier discussed in [TS], now all relative signs of the terms in (12.221) 
are identical, but the terms are otherwise unchanged! A proof of this transformation will 
be given in App. El 

To wrap up the above results, we would like to mention that the scaling functions 
fm(g), fo(g) and f(g) are part of a one-parameter family f K (g) interpolating between 
these three choices. The general function is obtained by multiplying the constants in 
( 12.191) by an overall factor of (1 + k). The resulting universal scaling function f K (g) is the 
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same as the sign-synchronized scaling function f(g) in (12.221) . but with the replacement 
C(2ri + 1) -> v^C(2n+l). 

We see that there are very interesting and seemingly natural ways to deform the 
scattering phase of [26, 5] while preserving Kotikov-Lipatov transcendentality. We will 
now argue in section [3] that the (conjectured) AdS/CFT correspondence, together with 
(conjectured) integrability and a (conjectured) crossing-symmetric strong-coupling phase 
factor indeed appears to single out one of the above choices. We will continue in Sec. 0] 
with the investigation of the kernels and scaling functions. There we will derive a closed 
form for the very same constants and a concise integral expression for the summed 
dressing kernel. 



3 An Analytic Continuation of Sorts 

In the following section, we shall investigate the constants 0} starting from string theory. 
For perturbative string theory it is useful to write the dressing phase 8kj in (I2.9|I2.10I) as 

oo oo 

9(u k ,Uj) = \] 2j c riS (g)(q r (u k ) q s (uj) - q s (u k ) q r (uj)) (3.1) 

r=2 s=r+l 

Here the excitation charges q r (u) are normalized as q r (u) = g r ~ 1 q r {u) differently from 
(12.101) . Consequently, the coefficient functions c rjS and (5 TjS are related by 

c r!S (g) = g 2 - r - s p r , s (g). (3.2) 

The strong-coupling expansion of c rs within string theory is non-trivial 



CrM = Y. c ^9 l ~ n - ( 3 - 3 ) 



n=0 



A proposal for the all-order strong-coupling expansion based on available data [19,32— 34J 
and crossing symmetry was made in Sec. 5.2 of 



r M (l-(-l) r+S )C(n) r[l( a + r + n-3)]r[l( a -r + n-l)] 

r ' s 2(-27r)T(n-l) 1 A ! r[|(a + r - n + 1)] T[\{s - r - n + 3)] ' 1 f 

This expression is formally 0/0 when setting n = 0, 1 so we have to use a proper regular- 
ization [23]: In these cases the correct expression is known from comparison to classical 
and one-loop string theory data [H? jl3"2Tf3"l] 

C W ('-'-"'1 (r -D(»-l) (35) 
n (s + r — 2){s — r) 

The expression for ci]} is easily recovered from (13. Ah through the limit ((n) / T(n — 1) — > 1 
when n — > 1. For n = we first set r, s to some integer values with s > r. If s > r + 1 
the only singular term is 1/ r(n — 1) which guarantees d°l = 0. This explains the 5 r+ i jS 
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term in (13.51) . For s = r + 1, however, the term r[|(s — r + n — 1)] diverges at n = 0. In 
this case we can rewrite the expression with arbitrary n as 

(») 2(n-l)C(n) £(l+|n) r(r - 1) T(r - 1 + \n) 

c r,r+i - i ~\ nTI ' , nr/o tv~ , i i„n ■ 



-27r)«r(n + 1)T(2 - hi) T{r + 1 



2 



n 



Here we can easily set n = and obtain cf]. +l = 1 as required by (13.51) . 

In order to compare to gauge theory, we have to find an expansion of Cr, s (g) at weak 
coupling according to (I2.10p . The main obstruction in the investigation of the function 
c rtS (g) is that the known expansion at strong coupling is merely asymptotic: For fixed 
r, s and even n the sequence of coefficients c^} terminates at n = s — r + 1 [23] . However, 
the odd-n sequence does not terminate and grows factorially as 

c (n+2) n 2 

' +0(n). (3.7) 



c. 



(n) 54^2 



r.s 



(n) 

Consequently, the series of c r J is asymptotic and has zero radius of convergence around 
g = oo. It is Borel summable, but the expression for the coefficients may be too complex 
to perform the sum in practice. 

In order to extrapolate to weak coupling, let us consider a simple model function 
first. In [23] some similarities between the phase and the digamma function *ff(z) = 
<9 2 logr(,2) were observed. This function has the following asymptotic expansion for 
large and positive z 

*(1 + z) = logz + ]T ^ , c n = -- S = (-1)»C(1 " n) , (3.8) 

n=l 

where B n are the Bernoulli numbers which can also be expressed through the zeta func- 
tion as B n = (— l) n+1 n((l — n). Conversely, the expansion around z = reads 

oo 

V(l + z) = - 7E + ^g fc z fc , c fc = -(-l) fc C(l + A:) (3.9) 

k=\ 

with 7e being Euler's constant. The curious observation is that the expansion coefficients 
for large z and for small z are almost the same. They are related by 

C n = -C- n - (3.10) 

In other words, we obtain the expansion coefficients by analytic continuation of the index 
parameter n to negative values. 

Could it be that a similar relation holds also for the expansion of the coefficients 
c r,s{g) a t strong and weak coupling, respectively? Admittedly, this is a very wild guess, 
but it turns out to be literally true. The expansion of c r>s (g) at weak coupling reads 

oo 

CrM = -Y. c ^ n) 9 1+n - (3- 11 ) 

71=1 
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In the following we shall not only provide indications for the mathematical correctness 
of this statement, but also argue that it leads to a consistent physical picture. 

As a first test, let us compute some c^} with negative value of n. Straight evaluation 
turns out to yield 0/0 in most cases and the expressions need to be regularized. We will 
therefore use the same procedure as in the cases n — 0, 1 and set r, s to their expected 
values first. Afterwards we take the analytic continuation of c^} for the expected value 
of n. In this way we obtain, e.g. 

4;3 1} =0, 4; 3 2) = -4C(3). (3.12) 

The absence of a n = — 1 contribution translates via £ — (r + s — n — l)/2to the absence of 
a £ = 5/2 contribution, where £ represents the number of loops in gauge theory. In general 
all contributions for odd negative n should be zero to guarantee absence of contributions 
at fractional orders. In that case the expansion coefficients at weak coupling would read 

m = -air- 2 *-*. (3.i3) 

The first non-trivial contribution appears al@ 

/3g = 4C(3). (3.14) 

This is encouraging for several reasons: Firstly, it influences anomalous dimensions start- 
ing from four loops. Currently, we know fljfl = for £ = 0, 1, 2 only, so this is perfectly 

(3) 

consistent with available data from gauge theory. Moreover, p\ 3 is the first allowed 
contribution according to (12.181) . In addition, the value 4£(3) has transcendentality 3 in 
the counting scheme of [H]. This is precisely the right transcendentality required for a 
uniform transcendentality of the scaling function [13] . Finally, the value (13.14p agrees 
precisely with the first in the series of sign reversing constants ( 12.21D ! 

Let us now rewrite the coefficients in a form more convenient for negative n. We use 
the identities 

C(l- z) =2(2ix)- z cos(\ixz)T(z)C ! (z) and T(l - z) = — -. , (3.15) 

2 sm(irz)T(z) 

and obtain for integer r, s a new form for the coefficients 

(n) (1 - (-l) r+s ) cos(f vm) (-I) 8 - 1 - 11 C(l - n) T(2 - n) T(l - n) (r - l)(s - 1) 



c r,t 



r[|(5 -n-r-s)} r[|(3 - n + r - s)\ T[±(3 -n-r + s)} r[±(l - n + r + s)] 

(3-16) 

The factor cos(|7m) makes it apparent that n must be even and that there are no 
fractional-loop contributions. Furthermore the factor 1/ r[|(5 — n — r — s)\ imposes the 



5 A curious observation is that the leading contribution to the function /?2,3 = C(3) A 3 /10247r 6 + . . . 
resembles the value 3C(3)A 3 /5127r 4 [3S] of a certain X-shaped diagram to a circular Maldacena- 
Wilson loop [37j (the conjectured value A 3 /4! 3 in [35] is almost true). If a contribution proportional to 
C(3) survives in the sum over all diagrams at this order, it would imply that the ladder approximation [38 
is not complete. This may seem unlikely, as this approximation yields a consistent strong-coupling 
result [55]. However, it is not yet excluded either that non-ladder contributions have a mild, but non- 
vanishing effect. 
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Figure 1: Left: Plot of h2s(z) with asymptotic value 8/37T. Right: Plot of the partial sums in 
c 23{g) with weak coupling expansion up to radius of convergence g = \. 

lower bound — n > r + s — 3 for the loop order £ = (r + s — n — l)/2>r + s — 2 consistent 
with (12.181) . As we discussed already in the previous section [21 

#2~C(2 + 2*-r-s) (3.17) 

this is precisely what is required for a uniform transcendentality of the scaling function. 

Being convinced of the physical plausibility of the proposal, we shall now investigate 
the equivalence of the weak and strong-coupling expansions for the choice r = 2, s = 3. 
The proposed weak-coupling expansion is 

oo 

c 2 , 3 (s) = -£4 2n y +2n (3- is ) 

71=1 

with 

_(-*0 = 4(-l)"C(2n + l) r(2n)T(2n-l) 
2 ' 3 " I» I> + 1) I> + 2) I> + 3) ' 1 ' ' 

We shall write £(1 + 2n) using its series representation 

l)»+i(2n-l)!(: 
(n-l)!n! (n + 1)! (n + 2)! 



OO OO 



M - V- V- 4 (- X ) n+1 ( 2n ~ ^ ( 2n ) ! (^/ fc ) 1+2n room 

n=l k=l 

The sum over n can now be evaluated explicitly to 

, , A, , , , , 1 l + 16z 2 , N 1-42 2 , . . . 

02,3(9) = Y, hiMk) , h 2>3 (z) = — + 3 K 4iz - — — E te , 3.21 

k=l 

where K and E are the elliptic integrals of the first and second kind, respectively^] The 
function fo 2 ,3, some partial sums as well as the series expansion of 02,3 are displayed in 

Fig. m 

It is now possible to connect to strong coupling. The expansion of the function /i 2 ,3 
at large positive z is 

h (z)= 8 1 I 31 ° g(lfo) 5 I 51 ° g(lfo) I 1 I f3 22) 

2,31 ; 3tt z Atcz 2 8nz 2 512tt2 4 2048ttz 4 " ' ' K ' } 



D We use the convention that the modulus appears in a squared form as (4zz) 2 
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Figure 2: Left: Plot of partial sums in 02,3 (<?) with asymptotic strong coupling expansion. Right: 
We have plotted xc2 t 3(l/x) instead to improve the view of the strong-coupling behavior. 



Zeta function regularizatiorfl of the sum yields 

Y^h 2 Mk) = g h 2 , 3 (l/z)dz+ 8uX>] U_J) 
Jo 



k=i 



3tt g 

3 C(-2) 5 C(-2) 5 CM) 1 C(~4) 
4vr g 2 8tt # 2 512tt c/ 4 2048vr c/ 4 
4 1 3C(3) 15C(5) 

5 3tt 12# 16vrV 2048vrV 

oo 

= E C &V-. (3.23) 

n=0 

Here the sum over fc n is regularized to ({—n) while the sum over k n log(fc) yields — C'( — n). 
We have verified the perfect agreement of the expansion with the strong-coupling coef- 
ficients C23 to order g~ 100 leaving little room for error. In App. O we present a proof 
of this statement to all orders. In the proof we observe some exponential corrections of 
the type 0(e~ 9 ) to the series. These corrections are not unexpected because the series 
is asymptotic. A plot of the partial sums and asymptotic strong-coupling expansion of 
c 2i 3 is given in Fig. [21 

We have also confirmed the same agreement for other values of (r, s): We can always 
replace the ((1 + 2n) in c r>s by a new sum and write c r>s in terms of a function h rjS (z) 
without odd-zeta terms 



fc=i 



7 Zeta function regularization is closely related to Euler-MacLaurin summation. However, the latter 
has to be adapted to deal correctly with terms k n log(fc) in the expansion of h, see [35] • We thank 
S. Frolov for providing us with this reference. 

8 This form of the phase may be suggestive of the presence of bound states [40] (or in general larger 
multiplets [41 j ) in intermediate channels during the scattering process. For instance, the dispersion 
relation as well as other characteristic quantities of a bound state of k particles is typically obtained by 
replacing g — > g/k in the expressions for a single particle. 
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The series in g can now be written formally as the hypergeometric function 



2cos[±vr(s - r - 1)1 T(r + s - 2) ^ 9 

KM = ^jr ( TvFT ^ " z ^ *; * 3 ' 25 

1 [r — 1) 1 (s — 1) 

with arguments a — [|(r+s — 2), i(r+s — 1), |(r+s — 1), |(r+s)] and 6 = [r, s, r + s — 1]. 
This hypergeometric function can be evaluated to the general form using elliptic integrals 

W») = + |iS K(to) + E(4fa). (3.26) 

Here P, Q, R are some polynomials of degree (r + s — 3)/2 with rational coefficients. The 
polynomials P and i? have special factors for all value of r, s, namely P(z 2 ) ~ ,2 2 ( r_2 ) 
and Q(z 2 ) ~ (1 + I62: 2 ). As in the above case of (r, s) = (2,3), we use zeta function 



regularization to evaluate the sums (I3.24f) at strong couplings We have further expanded 
some of the functions (13.25!) at strong coupling and they fully agreed with the proposal 
in [23] upon zeta function regularization in all cases tested (all pairs of r, s with s < 10 
to 20 orders). In App. Owe will finally show for general (r, s) that c r ^ s (g) at large g 
yields the first two orders cl°} and cl]} in (13.51) which are known from perturbative string 
theory [19j[33]. 

It is therefore fair to claim that the strong-coupling expansion (I3.3P and the weak- 
coupling expansion (13. lip describe one and the same function c TyS (g) for all (r, s). How- 
ever, we should stress that the weak-coupling expansion is more versatile than the strong- 
coupling one: The functions h rjS have singularities at g = ±4, and thus the coefficients 
c r . tS (g) have singularities at g = ±jk for all positive integers /c0 Consequently, c TjS (g) 
has a finite radius of convergence around g = 0, and the weak-coupling series defines the 
function unambiguously. In contradistinction the strong-coupling expansion is merely 
asymptotic and does not fully define the function. 



4 Magic Kernels 

Now we will return to our study of the integral equation for the fluctuation density in 
the presence of a perturbative dressing phase 9 7^ 0. As noted in section [2} there appears 
to exist a very special choice for the constants /3^fj in (I2.10p such that the odd-zeta 
contributions to the scaling function based on a trivial phase 9 = (f m (g), see (12.81) ) 
either cancel {fo{g), see (12.201) ). or synchronize their relative sign w.r.t. the even-zeta 
contributions (/((?), see (I2.22p ). The constants of the latter case appeared to differ from 
the ones of the former case by a factor of two. Let us now prove these statements, and 
derive the relevant set of constants (3^% 

An observation of central importance for this section is that the simplified scaling 
function fo(g) may be obtained from an effective kernel much simpler than (I2.13|2.14p . 

9 The polynomial P generates even loop orders in string theory, while Q, R generate the odd ones. 
Its coefficients are thus directly related to the coefficients c^f 1 . 

10 The singularities arise from the elliptic integrals in (|3.26p and therefore appear to be related to a 
degeneration of the complex tori defining the kinematical space for the bound states in [301I3T] . 
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One just needs to decompose the "main scattering" kernel ( 12. 6ft into two parts 

k m (t,t') = k (t,t') + k 1 (t,t'), (4.i) 

which are even and odd functions, respectively, under both t — > —t, t' — > —t'. Explicitly 
these functions read 

{i a t m Jot?) - 1> j (t) m?) 

K {t, t ) = — — , 

f> ( . t' .hit). J {t')- t. J {t). hit') 

K 1 {t,t) = _ — . (4.2) 

It turns out that the component K\ is responsible for all odd-zeta contributions to / m . 
If we eliminate it and set K = K in (12.21) we obtain just the same functions o"o and 
fo(g), see ( I2.20p . Note that again ( 12. 4ft remains true under this symmetrization, and 
( 127711 holds with K m -> K . 

We can use this observation to obtain an almost closed expression for the dressing 
kernel. We start from the integral equation with a purely even kernel (14. 2 j) 

a (t) = ^— AT (2pt J 0)-4p 2 / dt' K (2gt,2gt')a (t f ) , (4.3) 

which leads to fo(g). Replacing K by K m — Ki, using Ki(t, 0) = 0, and substituting fr 
once into the ensuing second convolution only, we derive the relation 

<r (t) = —^-\k m (2gt,0) + k c (2gt,0) 



e* - 1 



with the function JC. 



-4/^ dt' (k m (2gt, 2 gt') + K c {2 gt,2 gt')^ a {t')\. (4.4) 



K c (t,t')=4g 2 dfR^gt'^^—K^gtht'). (4.5) 
Jo e — l 

Now (14.4p is just the original integral equation (12.21) . but with the kernel K = K m +k c . 
Therefore we can interpret K c as a dressing kernel to generate the simplified dressing 
function fo(g). In fact, it agrees precisely with the expansion of the general dressing 
kernel (12.13112.141) using the first few coefficients found in (I2.19p ! The expression for the 
dressing kernel K c allows us to derive a closed expression for the constants (3^) leading 
to cancellation of odd-zeta contributions. The details of this derivation are presented in 
App. 13] the final result is 

Surprisingly, these coefficients agree, up to a global factor of |, with the "analytic 
continuation" of the conjectured strong-coupling dressing phase constants (I3.13|3.16p ! 
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-1 -0.5 0.5 1 

Figure 3: Top to bottom curve: Plot of the scaling functions f m (g)i fo(d) an d f(g) with 
9 = lV%- 

The results of the previous section then suggest that the correct crossing-invariant choice 
of constants is just twice (}4.6j) . Of course, these constants agree with those in (I2.2ip . 
The corresponding kernel for the integral equation (12.21) of the universal scaling function 

f( g ) is, cf. (izara 

k(t,t f )=K m (t,t , ) + 2K c (t,t'). (4.7) 

What does this imply for the universal scaling function? We "experimentally" ob- 
served in section [21 to rather high order, that the doubled constants (14.81) lead to sign- 
reversals in those odd-zeta terms which do not have the same relative sign as the even-zeta 
terms in the original proposal (12. 8p . Up to this minor modification, the associated scaling 
function agrees with the one based on a trivial phase. We will present a proof of this 
statement in App. [B] 

It is interesting to investigate the analytic structure of the three candidate scaling 
functions. In each case one would like to find the locations of all singularities in the 
complex (yf-plane, and determine wether they correspond to poles, branch points, or are 
of essential type. The convergence properties of the weak coupling expansion series are, 
in particular, clearly related to the singularity with smallest \g\. 

The integral form of the dressing kernel K c allows to derive conveniently the first few 
orders in the expansion of the scaling function. The analytic expression is given (I2.22p . 
Beyond 20 orders, however, the analytic calculation slowed down substantially due to 
an exploding number of partitions for the odd-zeta contributions. This complexity can 
be reduced by reverting to a numerical computation. We were able to produce all three 
scaling functions f m (g), fo(g) and f(g) up to 50 loops, O(g wo ), at 1000-digit accuracy 
(needed for intermediate steps). A plot of these three functions is shown in Fig. [31 

We have then used the quotient criterion to investigate the convergence properties 
of the series. Our results indicate that the radius of convergence in g is ~ in all three 
cases. The related singularity is situated at g — ±| or, put differently, at A = — n 2 . 
The accuracy of this result was very high and the position matches with the singularity 
of the function in (I3.2ip . The exponents of the singularity / ~ (A + 7r 2 ) a appear to be 
a = +1, +|, +| for the three functions f m (g), fo(g), fig), respectively. We were not 
able to produce a very accurate result in this case, we expect it to be accurate within a 
few percent. 
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In conclusion, the final proposal for the, hopefully correct, expansion coefficients of 
the dressing phase (12.10112.1 ip reads 



The resulting universal scaling function f(g) in a weak-coupling expansion was presented 
in (I2.22p to the first few orders. Its convergence properties seem to indicate a singularity 
at A = — 7r 2 with exponent +|. 

The strong-coupling expansion of the scaling function f(g) remains to be evaluated 
within the Bethe ansatz framework. A comparison to available data [4*21, |4"3] from string 
theory would constitute an important test of our proposed dressing phase. The leading 
order at strong coupling based on the AFS phase "T"*"] (which is consistent with our 
proposal) indeed agrees jS] with an explicit string theory calculation [12]. Conversely, 
in the case of the un-dressed function f m (g), the integral equation is troublesome and 
leads to a fluctuating solution [29]. Hopefully the situation is improved by the correct 
dressing phase in order to find agreement with string theory. 



5 Conclusions and Outlook 

In this paper we have derived an expression (12. 10112. 1 1|I4.8[) for the complete weak-coupling 
expansion of the dressing phase 9 in the Bethe ansatz [H] of planar M = 4 gauge theory. 
It is based on a crossing-symmetric proposal in Sec. 5.2 of [23] for the strong-coupling 
expansion of the dressing phase in conjunction with some curious structures of cancella- 
tions and replacements when a dressing phase modifies the scaling function f m (g) based 
on a trivial phase, as originally derived in [13]. Whether or not our proposal describes the 
AdS/CFT system accurately, we have shown that there exists a natural expression for 
the dressing phase which appears to interpolate between currently known spectral data 
at weak and strong coupling. We are hopeful, but have not shown, that any alteration 
of our proposal is likely to violate the transcendentality principle in gauge theory and/or 
crossing symmetry in string theory as well as leads to potential contradictions with the 
available data. In that sense our proposal lends support to an exact interpretation of 
the AdS/CFT correspondence [17], at least in the case of M = 4 gauge theory in the 
't Hooft limit and free IIB superstrings on AdS^ x S 5 . 

The dressing phase 9 is given as a multiple expansion in excitation charges and the 
coupling constant. The radius of convergence w.r.t. the 't Hooft coupling A around 
A = is |A| < 7T 2 . The singularities are situated at A = —n 2r K 2 for all non-zero integers 
n. Consequently, perturbative gauge theory appears to unambiguously define the phase 
and its analytic continuation. Note however that we have not succeeded in summing over 
the modes r, s of (3 r ,s{g) or even study the analytic structure of the phase 9 in general. 
Conversely, the strong-coupling expansion of the phase can merely be asymptotic because 
A = oo is an accumulation point of singularities. The asymptotic string expansion 
proposed in Sec. 5.2 of [23j seems to be in perfect agreement with our proposal, but it is 
not sufficient to uniquely define a complex function. 
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Two methods of testing our coefficients come to mind: The first is a direct computa- 
tion of the scaling function f(g) at four loops. Until recently, it was known up to three 
loops [niEj, see also [H], but with methods based on the unitarity of scattering amplitudes 
it is within computational reach [H1H2]- Our prediction (3^1 = 4/3 = 4£(3) leads to a 
four-loop contribution to the universal scaling function 

/(*) = . ..-16 Q^ 6 + 4 C(3) 2 ) 9 8 +... 

«... - 3.01502 x 1(T 6 A 4 + ... . (5.1) 

The final outcome of the four-loop calculation by Bern, Czakon, Dixon, Kosower and 
Smirnov [TJ], which was performed in parallel and completely independently from our 
analysis, is 

f( g ) = ... - 64 x (29.335 ± 0.052) g 8 + . . . 

= (3.0192 ± 0.0054) x 10~ 6 A 4 + . . . . (5.2) 

The perfect agreement between the two numbers provides further confidence that our 
proposed dressing phase correctly describes the asymptotic spectrum of anomalous di- 
mensions in planar M = 4 gauge theory. At the least it ascertains that the dressing 
phase is non-trivial starting from 0(g 6 ), and that the constant for the first correction 
in [13] takes the value (3 = = C(3)- Note, however, that we have to assume the 

exact universality of the scaling function, i.e. that it is exactly the same for twist-two 
(as assumed in [12J) and for higher-twist (as required for the asymptotic Bethe ansatz). 

A computation of the function f(g) to even higher loop orders would provide more 
stringent tests of our proposal. Given the remarkable success of unitarity-based methods 
in [6ll7lfT2] we do not dare calling this impracticable. However, some yet-to-be-discovered 
iterative structure is likely to exist if our comparably simple integral equation holds true. 
With many perturbative orders available an extrapolation to strong coupling may be at- 
tempted and compared to string theory results [^2"jH3"] via the AdS/CFT correspondence. 
Such an analysis was indeed performed in [12] using the KLV |45j and Pade approxi- 
mations. It is based on their independent guess of the function f(g) in (12.221) and how 
to obtain it from the earlier proposal (12.81) by synchronizing all signs. Strangely, their 
approximations, based on using the first seven or eight terms for example, appear to 
indicate that (12.221) is about 5% away from the expected strong-coupling value of [42J. 
Non-coincidence of the values would be rather puzzling with regard to the discussion 
at the end of Sec. H] and Sec. [3l Would it imply that f(g) at strong coupling does not 
describe the string states in [42}H3]? Are there some non-perturbative corrections? To 
that end, it would be worth investigating whether an adapted scheme which takes into 
account all singularities at A = —n 2 7i 2 and their accumulation at A = oo can produce an 
extrapolation of f(g) closer to the expected string theory values [12HS]. 

Another way of testing the proposal involves local operators of small classical di- 
mension. A gauge Bethe ansatz modified by a dressing phase changes the anomalous 
dimensions A of all operators in all sectors. E.g. in the su(2) sector we would find for 
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the length L = 5 operator TrX 2 Z 3 + . . . (this case is actually equivalent to the sl(2) 
twist=length three operator Tr D 2 Z 3 + . . .) to four loops 

A = 5 + 8/ -24s 4 + 136s 6 - 16 + 2^ 3 3 } ) g 8 + . . . . (5.3) 

For the specific case of constants (14.81) we then have 

A = 5 + 8 g 2 - 24 g 4 + 136 g 6 - 16 + 8 C(3)) /+.... (5.4) 

Intriguingly, this predicts the appearance of a non-rational term in the four-loop anoma- 
lous dimension of a finite length operator! Such contributions are impossible if the 
dressing phase is trivial, i.e. if 9 = 0. They are certainly not excluded, and a priori even 
likely, from the point of view of perturbative field theory. This would then also rule 
out the all- loop "BDS" ansatz originally proposed in [26J, and as a consequence also its 
finite-length description through the Lieb-Wu equations of the Hubbard model, cf. (68) 
in [IS]. It also eliminates the Hubbard Hamiltonian as an exact candidate for the su(2) 
dilatation operator [46J. However, a gauge theory dressing phase (I2.10p whose constants 
preserve transcendentality for twist operators, such as in the specific suggestion (14. 8ft . 
would have the following intriguing property affecting operators of finite length and spin: 
Up to (at least) wrapping order, the BDS ansatz would give the correct "rational part" 
of all anomalous dimension matrices, and the Hubbard Hamiltonian would emulate the 
"rational part" of the su(2) dilatation operator to (at least) wrapping order. 

An interesting question concerns the computation of anomalous dimensions beyond 
wrapping order. The simplest case is the four-loop anomalous dimension of the Konishi 
field, i.e. the su(2) length L = 4 operator TrX 2 Z 2 + . . ., or, equivalently, in sl(2) the 
lowest spin S — 2, length L = 2 operator Tr D 2 Z 2 + . . .. If we, a priori "illegally", apply 
the asymptotic Bethe ansatz ( 12.9112.101) we find the four-loop result 

A = 4 + 12 g 2 - 48 g 4 + 336 g 6 - 16 ( 7 -f + § g 8 + . . . . (5.5) 

(3) 

Recall that (3^ = is the "BDS" case with dressing phase 6 = 0. With our conjecture 
(14. 81) we then have instead 

A = 4 + 12 g 2 - 48 g 4 + 336 g 6 - 16 (?f + 18 C(3)) /+.••• (5.6) 

However, as opposed to (15.41) . this result is much less certain. It would only be true if our 
asymptotic Bethe equations turn out to be exact. Conversely, recall that the Hubbard 
model [10] yields for this state 

A = 4 + 12 g 2 - 48 g 4 + 336 g 6 - 16 x 318 g 8 + . . . . (5.7) 

Could it be that the Hubbard model still properly yields the "rational part" of fields 
even beyond wrapping order? This might then result in 

A = 4 + 12 g 2 - 48 g 4 + 336 g 6 - 16 (318 + 18 C(3)) /+..., (5.8) 

but, in the absence of a yet-to-be constructed rigorous Bethe ansatz for the finite system, 
this is of course a quite unfounded speculation. 



17 



Of course, there are more general questions: It would be very valuable to find a 
closed form for the proposed dressing phase 9 and to prove its correctness. Likewise, the 
exact integrability of the planar AdS/CFT model still lacks a rigorous proof. Are the 
Bethe equations at finite coupling and finite length exact or is there a more fundamental 
description of the planar spectrum of which the Bethe equations are merely some limits? 
For instance, the exactness may be flawed by certain terms which are exponentially sup- 
pressed at long states [UJ . The "covariant" approaches started in P8 | H6 | H9] may provide 
an answer to this question as well as recover the dressing phase from a more fundamental 
(and simpler) description. Other promising approaches are the thermodynamical Bethe 
ansatz (see [50]) and Baxter equations (see [5"T]). 

It is curious to see that in our proposal the deviations from a trivial phase start to 
contribute to anomalous dimensions at four loops, which is a rather high perturbative 
order. If true, this might indicate that intuition gained from some perturbative com- 
putations, even if they extend to two or three loops, may be deceiving. For instance, 
a non-zero value for (3^1 = 4/3 leads to a four-loop breakdown of weak-coupling BMN- 
scaling [2] in Af = 4 gauge theory. This possibility was discussed by several authors in 
the past in order to reconcile discrepancies between gauge and string theory structures 
at, respectively weak and strong coupling. See in particular [52j[T6j[26j[53l[32l[5l]. Our 
proposal suggests that this is indeed the case. The assumption of [2] that the BMN limit 
leads to a dilute gas with weakly interacting particles is not compatible with our phase. 
It should be noted that the perturbative dispersion relation obtained in [55] is indeed 
compatible with our proposal. However, it cannot be translated directly into anomalous 
dimensions of local operators as suggested in [21155]. 

The same applies for the spinning strings proposal [56J : The weak-coupling expansion 
in the effective coupling constant A' ~ A/ J 2 leads to terms of the type (A') 4 J 2 which 
diverge in the limit of large spin J. The investigation of spinning string solutions from 
the gauge theory side started in [57] will therefore fatally break down at the four-loop 
level. Nevertheless, the proposals [2,56j have been of tremendous importance to initiate 
the study of stringy and integrable structures in J\f = 4 gauge theory. 

Finally, it is worth pointing out that ^2,3(9) — @(9 6 ) is reminiscent of the findings [SH] 
within the plane wave (BMN) matrix model [2]. Also there, a non-trivial dressing phase 
was observed starting at the same order. The main difference is that the coefficient {3\ 3 
for the BMN matrix model is rational while ours is a rational multiple of C(3). This may 
potentially be related to the fact that the matrix model has only finitely many degrees 
of freedom while M = 4 SYM is non-compact. 
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A Expansion of the Magic Kernel 

In this appendix we present a weak-coupling expansion of the magic kernel K c in (14.51) 
leading to an expression for the coefficients (5r) s of the dressing phase. 

A useful integral representation for the even and odd parts was found in [IB] : 



K (t, t') = [ dX A J (A t) J (A t') , 
Jo 

k 1 {t,if)= [ dxx Ji(xt) Ji{\t') . (a.i; 

Jo 



'0 

With their help, we may calculate K c to be 



/•OO /•! j II 

K c (t,t') = 4g 2 dt" / dX 1 X 1 J 1 {X 1 t)J 1 {X 1 2gt")— l - 

Jo Jo e — 1 

rfA 2 A 2 J (A 2 2(7t") J (A 2 f) 



o 



oo 



-4^(-l)V M+1 C(2/i+l) (2//)! x 

E •■• - — ^ 



ra=0 



n\ (n + 1)! (/i — 1 — n)! (/i — 1 — n) 



[ 1 dX 1 X 2 1 {n+1) J 1 (X 1 t) f 1 dX 2 Xl i ^ n) - 1 J (X 2 t') (A.2) 



where we have first expanded the Bessel functions with argument Ai j2 2 g t" in g and then 
executed the t" integration. The two remaining parameter integrals in the last line may 
be found from a standard formula. In the first integral the power of Ai in front of Ji(Xit) 
is always even, and in the second one the power of A 2 in front of Jo(X 2 t') is always odd. 
Both integrals therefore yield finite sums of Bessel functions: 



f dXxX 2 ™ JxiXxt) 
Jo 

(m " 1)!m ' B-i)< (^\)^« 

fc=i ^ ' 



(2m) 



f dX 2 Xl m - 1 J (X 2 t') 
Jo 

(m-l)!(m-l)! " (2m - 1\ 21 - 1 
(2m -1)! V m - / ) ~P~ } " (A ' 3) 



On substituting these into ( 1A.2I) we can cancel the factors which depend on n but not 
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on k, I. We re-order the n, k, I sums to find: 



oo 



K c {t, = V M+1 C(2/i + 1) (2//)! x 

/x=l 

j£+K/l+l 

£ (-l) fe +^(2fc) J 2fc (i) (2/ - 1) Ja.!^) x 
fc,z>i 

n 5i + 1 - + 1 + - ^ - + ^ - 1 - ' ( ' A ' 4 ' ) 

The sum over n can be done in closed form: 

/it— i 



n 5i { n + i - k y-( n + i + k v-{^- i - n V-(^+ 1 - 1 - n V- 



(2/x + 1)! V/" + 1 - fc - V VA 4 + 1 + k ~ 1 
Collecting results, we get 

k+l<n+X 

tt 

{2k) {21 -I) ( 2/1+1 \ ( 2/i + l 



A 

(U=l fc,Z>l 



2// + 1 V/i + l- fc-ZyV/i + ! + £;-/ 



C(2//+l). (A.6) 



This is precisely of the form of a kernel stemming from a non-trivial dressing phase 9, 
as we can see by comparing to (I2.14p ! We can even read off the constants f3^}: 

Finally, instead of using the integral representation flA.ip . we may also cast the two 
kernels Kq, K\ in the form of infinite sums0 

2 °° 

Ko{t, t') = — ^(2n + 1) J 2n +i{t) J 2 n+i(t') , 



n=0 

oo 



KiM = ^Y,^n) J 2n {t) J 2n {t') . 



tt' 

n=l 



It follows that 

*c(«, = ^ ]T(-l) fc+ W, 2 / +2 (s) J»(f) J a+1 (0 (A.8) 



2 ' ' 



fe=l /=o 



11 This result was first obtained by D. Serban 
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with the coefficient functions 



CrM = 2 C os[^(s - r - 1)] (r - - 1) dt J ^t) J s -^(2 g t) ^ 

We can use this form to present an alternative proof. We firstly expand l/(e* — 1) into a 
geometric series and secondly rescale the integration variable t by g such that the integral 
will manifestly depend on the ratio g/n only 

oo 

Cr,s{g) = ^hr^g/n), (A.10) 

n=l 

f°° dt 

h r>s {g/n) = 2 cos[±vr( S - r - 1)] (r - l)(s - 1) / - exp(-tn/g) J r . x {2t) J,_i(2t) . 

Jo t 

The integral can now be performed and yields precisely the hypergeometric function 
derived earlier in ( 13.25(1 . 

The discussion of the strong-coupling behavior of the dressing phase is tantamount to 
an analysis of the strong-coupling limit of the integrals c rjS (g). In App. Owe present, by 
way of example, the full expansion in the case (r, s) = (2, 3). We also show that (IA.10I) 
reproduces the AFS phase [19] as well as the Hernandez-Lopez correction [33] . 



B Flipping Odd-Zetas Contributions 

Here we show that the replacement of the kernel K m — > K m + 2K C leads to the simple 
replacement ((2n + 1) — > i((2n + 1) in the scaling function f m (g) — > /(<?)• 
Let 

*(*) = -* (B.l) 
e { — 1 

For any kernel K(2gt,2gt') and any function f(t) we define the product 

POO 

(K * /) (t) = / dt' K{2 gt,2g f) z{?) f(t') . (B.2) 
Jo 

The integral equation may be written as 

a(t)=z(t)(^k + k c )(2gt,0)-4g 2 (^K + & + K c ) *~) (t)J . (B.3) 

In the following we need not consider the internal (/-dependence of the kernels. Within 
the radius of convergence of the iteration the solution of any Fredholm equation of the 
second type is the resolvent 

a(t) = z(t) jr (-4g 2 T {[{Ko + K 1+ K c ) *]» p\ (t) . (B.4) 

n=0 V ' ' 
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Here P{t) = (K +K c )(2gt,0) and K c (2gt, 2gf) = Sg 2 ^ * K )(2gt, 2gf). The Fredholm 
resolvent must have the form 

a(t) = z(t)J2(-±9 2 TVn(WnKa)(t) (B.5) 

n 

because the dressing kernel K c consists itself of Kq, K\ which are convoluted by the star 
operation (and it also carries the appropriate extra power of the coupling constant). By 
n we denote a string built of 0, 1 with total length n, and Wn is the corresponding string 
of K Q , K\ acting on the potential Kq by the * product defined above in (1B.2H . The ry„ are 
numerical coefficients. The integral equation ( IB. 31) determines these numbers uniquely: 
We substitute (1B.5I) into the equation and pick the 0(g 2( - n+2 ^) terms. All words are 
assumed independent; in particular the four double iterations of an arbitrary word Wn- 

WW WW W 1>0 ,n, W 1X n- (B.6) 

This yields the four equations 

V0,0,n = — V0,n , 
V0,l,n = ~ r ll,n j 
Vl,0,n = -V0,n ~ 2Vn, 
r]l,l,n = ~Vl,n ■ 

The first two equations may be summarized as 

V0,rn = -Vm (B.7) 

from which it follows that the third equation is: 

Vl,0,n = +V0,n- (B.8) 

Recursively, these equations imply that all the rj coefficients are merely signs. Second, 
( IB. 81) implies that the operation of appending K\ after Kq is the only one that induces 



a sign flip as opposed to the iteration of the main scattering kernel K m = Kq + K\ alone 
(in that case the second term in the r.h.s. of (1B.7[) would be absent so that 771^ = —rjm, 
too). 

Next, let 

£ = {t 21 : I E No}, O = {t 2l+1 : / E N }, (B.9) 

denote the even and odd powers of t, respectively. In order to obtain the scaling functions 
at weak coupling we will eventually expand the kernels K 0t i(2gt,2gt') in the coupling 
constant; Kq (Ki) is even (odd) in both arguments, respectively. Let us consider any 
one term in the expansion of a kernel convoluted on a power of t by the star product. 
We observe the graded structure 

C(2m) t 2l \ 
((2m + l)t 21 ' 
C(2m + l)t 
C(2m) t 2V+1 



Kq : 


£ - 


^£ : 







O - 


^£ : 


0+1 




£ - 


: 


t 2l 




O - 


: 


f 2l+l 



(B.10) 
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The first integral in a chain of convolutions always acts on an element of 8 because the 
potential is even. Likewise, the energy integral adds a final convolution on the even 
kernel K (0,2gt). The contribution of a word W R to the scaling function thus contains 
one odd zeta function for each beginning and each end of a K\ string in Kq * Wn* Kq. In 
particular, if there are m such strings, the word gives a contribution with 2m odd zeta 
functions. 

Finally, the equations (IB.7IIB.8I) say that each creation of a K\ string induces a sign 
change as opposed to the iteration of the main scattering kernel on its own. Our result 
is: The scaling functions for both cases contain terms with 2m odd zeta functions. The 
relative sign of such terms is (— l) m . 

Note that the absolute sign of the terms with odd zeta functions cannot easily be 
fixed without explicit computation. Words with n insertions of K\ pick up a factor 
( — l) n from the g 2 expansion relative to those with K alone. Presumably, words with 
the minimal number of K\ kernels give the dominating contribution: K\ starts at 0(g 2 ) 
thereby restricting the number of ways one may choose the powers in the various kernels 
in a chain, that is reducing the number of similar terms. If this is so, the iteration of 
the main scattering kernel should indeed produce a factor of (— l) m for a term with 2m 
odd zeta functions relative to the leading purely even part. Consequently, the dressing 
kernel ought to align all signs. 

C Strong-Coupling Expansion 

In this appendix we give an exact derivation of the strong coupling expansion of 02,3 (g). 
We modify a technique developed in [20] for the discussion of the ground state energy of 
the half-filled Hubbard model. 
According to equation ( 1A.10I) . 



with 



f'(z) = Az-Az 2 F!(-i |; 2; z~ 2 ) - — 2 F X (A §; 3; z~ 2 ) 



1 

2' 2' ~ > 2^ zii ^2'2 ; 

/(z) = 2z 2 (l- 2 F 1 (-|,i;2;^ 2 )). ' (C.2) 



We want to use the residue theorem to express the sum as a contour integral. For any 
holomorphic function f(z) 



in 



h 2 (A7igz) {^g) 2 V 
because the terms in which the derivative does not fall on f(z) cancel out@ The integral 



r6Sm/49 sinh 2 (47r^) " (ing) 2 f [49) (C ' 3) 



12 The original article [60) uses a simple sinh in the denominator. This would lead to an alternating 
sum of residues, though. 
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Figure 4: Contour C and its decomposition (|C.5p . 



is thus written as 

° 2M ~ J c sinh 2 (4^) " I smh^gz) (1 ~ 2Fl( ~ 2 ' 2 ' 2 ' * )] ' {CA) 

The region encircled by the contour C contains all the points in/Ag up to n < Ag/e and 
f(z) must be holomorphic on it. Let e > infinitesimal. We shall consider integration 
over the closed contour, see Fig. 0] 

c = d u c 2 u c 3 u c 4 u c 5 , 

C x = (l/e)e i[0 ^ , 

C 2 = [-1/e, -1 - e] U [+1 + e, +l/e] , 
C 3 = (-1 + ee^M) u ( +1 + ee i[o !7 r]) ? 

C 4 = [-1 + e, -e] U [+e, +1 - e] , 

C 5 = ee* 10 '" 1 , (C.5) 

Note that the hypergeometric function in f(z) has singularities at z — ±1 and we shall 
deform the standard cut on the interval [—1,1] such that it lies below the contour C4. 
Let us first consider the integral over the large semicircle C\. Here we have 

/(*) = 1 + 0(1/*) (C.6) 

and for large enough 1/e we find (neglecting corrections in e) 

2 f dzf(z) 2 f dz _ _ (o) 



8719 / • J (A \ = / — 7YT. r = 9 coth^/e) = g = g. (C.7) 

J Cl smh {Airgz) J Cl smh {Airgz) 

The first term in f(z) leads to no singularity on the real axis and can be integrated 
straightforwardly over the remaining contours 

lim / l^9 2 dzz> = r° I6ng*dzz* = J_ = eg _ 

m -*°° J c 2 uc 3 uc 4 uc 5 sinh 2 (4vr5(z) 7„ 00 sinh 2 (4vr5(z) 125-5- 

For very small z the remaining hypergeometric function goes as 

4? 

^H,!^;,- 2 )^— (C.9) 
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so that the semi-circle C5 at the origin contributes 

16ng 2 dzz 2 11.^-2, 4 _ JD (r im 

c, sinh 2 (4^) * F ^-2>-*> 2 > z )-"^- C ^- ( C - 10 ) 

Since the integrand is an even function of z, the other three pieces add up to the real 
part of an integral along the positive semi-axis only: 

H{ 9 ) = -16V / . .f,f ; 2 F 1 (-|, §; 2; z' 2 ) 

Jc 2 uc 3 uc 4 smh (4tt5(z) 

= ^'*f i^fe'^-*^) (Cll) 

The real part of the hypergeometric function is easily read off from its one-parameter 
representation 

4 /•min(z,l) 

Re 2 F 1 (-i,i;2;^- 2 ) = - / dt VT^Vl - t 2 Z - 2 . (C.12) 

71 Jo 

In if (g) we readily estimate the contribution from z > 1 to be of order (D(e~ g ) due to 
the suppression by 1/ sinh(47rg,2) 2 . The remainder is 

H{g) = -128g 2 4r 7/ dt VT^Vl - t 2 z~ 2 + C(e" 9 ) . (C.13) 

Jo sinh (Angz) Jo 

To deal with these integrals we successively 

• expand y/l — t 2 into a series in t 

v m=0 v 



• calculate the t integral in each term 



f dt t^VT^^ = ^ 3^1) 

7o 4 r(m + 2) 



(C.15) 



evaluate the z integral. Here we may extend the range to the whole positive semi- 
axis since the error is again seen to be 0(e~ 9 ) 



I 



1 dzz 2m+3 ,00 dzz 2m +3 = 4 T(2m + 4)C(2m + 3) _ 

o sinh 2 (47n^) J sinh 2 (47n^) 16 ' ~ ~ (8ng) 2m + 4 [C h 

(C.16) 
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On collecting terms: 

R(fA = V r(m-i)r(m+|)r(2m + 4) C(2m + 3) 
^ r(m+l)r(m + 2) 7r 2 (87r^) 2m + 2 1 J 



m=0 

oo 



v 2 CH r[i(n + 2)]r[in] 

^ (-27r)»r(n - 1) r[i(6 - n)] r[j(4 - n)] * 



oo 



^c&V^ + OCe-*). (C17) 



n=3 



The second line is obtained by using some identities of gamma functions and setting 
n = 2m + 3. In conclusion: 

(2) oo 

C2M =c^9 + eg + °-f + H(g) = + O(e^). (C.18) 

We proved the correctness of the strong-coupling expansion coefficients c$ from equa- 
tion (13.4)) (see also (I3.23P ) for the case (r,s) = (2,3). The role of the zeta function 
regularisation used in the main text is simply to dispose of the terms of order 0(e~ 9 ) 
that become visible in this proof. 

The argument presented in this appendix should remain applicable for generic values 
of (r, s) . We hope to return to this general case in future work. We can however perform 
the analysis for the first two orders. By Euler-MacLaurin summation (or zeta function 
regularization), the leading two orders of c r}S (g) takes the form 



JO) 



h r , s (l/z)dz, cH=C(0)M«0- ( C - 19 ) 



In fact, one can show that the above argument yields the same answers. The integral of 
h r>s in the form (13.251) can be performed in closed form for all (r, s) 

JQ) 8(r-l)(s-l)sin[7r(r-s)] 

r ' s 7r ( s -r-l)(s-r + l)(s + r-3)(s + r-l) r+1,a 1 ' 

and one obtains full agreement with the AFS phase (13.51) [19] using that r, s are integers 
with r > s + 1. For the next order one can use the representation ( 1A.10I) of h r ^ s (z) at 
z = oo 

f°° dt 

eg =C(0)2cos[|7r(s-r-l)](r-l)(s-l) / - J r _i(2t) J._i(2t) 

Jo 1 

2 cos[|7r(s — r — 1)] 



cos[|7r(s-r-l)](r-l)(s-l) 2 
l_(_l)r+- ( r — — 1) 



2 J ' " ' 7r(s — r)(s + r — 2) 



7T (s — r)(s + r — 2) 
This is again in full agreement with the HL phase (13. 5p [33J. 



(C.21) 
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